Stochastic dynamics of a sheared granular 

medium 



Alberto Petri^ Andrea Baldassarri^ Fergal Dalton^ 
Giorgio Pontuale^ Luciano Pietronero^'^ 
Stefano Zapperi^ 

October 16, 2008 

^ CNR-Istituto dei Sistemi Complessi, Via del Fosso del Cavaliere, 00133 
Roma, Italy 

^ Dipartimento di Fisica, Sapienza, Universita di Roma, P.le A. Moro 2, 00185 
Roma, Italy 

^ CNR-INFM, SMC, Dipartimento di Fisica, Sapienza Universita di Roma, P.le 
A. Moro 2, 00185 Roma, Italy 

Abstract 

We experimentally investigate the response of a sheared granular medium 
in a Couette geometry. The apparatus exhibits the expected stick-slip mo- 
tion and we probe it in the very intermittent regime resulting from low 
driving. Statistical analysis of the dynamic fluctuations reveals notable 
regularities. We observe a possible stability property for the torque dis- 
tribution, reminiscent of the stability of Gaussian independent variables. 
In this case, however, the variables are correlated and the distribution is 
skewed. Moreover, the whole dynamical intermittent regime can be de- 
scribed with a simple stochastic model, finding good quantitative agree- 
ment with the experimental data. Interestingly, a similar model has been 
previously introduced in the study of magnetic domain wall motion, a 
source of Barkhausen noise. Our study suggests interesting connections 
between different complex phenomena and reveals some unexpected fea- 
tures that remain to be explained. 

1 Introduction 

The response of a granular material to shear stress has been widely investigated 
in the past. The first relevant report on the subject was perhaps that of Bag- 
nold [T] , pointing out a transition from sliding to stick-slip motion. Subsequent 
work (see e.g. references in [2]), was limited essentially to situations with more 
continuous or constant shear rate. However, we may observe that the response 
of a granular medium to an applied shear stress is often non-continuous, but 



takes place in the form of intermittent, erratic events, particularly when the 
applied shear rate is below some threshold. Despite the fact that this situation 
is frequently encountered in nature (earthquakes are a large-scale example), sys- 
tematic investigation of the statistics of fluctuations seems limited. One recent 
step in this direction was taken by Nasuno and coworkers U who analyzed the 
stick-slip behavior of a sheared granular medium, but considering a highly pe- 
riodic motion regime. Molecular dynamics simulations of a similar setup had 
been previously performed by Thompson and Grest [5]. The aperiodic regime 
was investigated by Dalton and Corcoran in a slightly different apparatus [31 [5] , 
who presented the statistics of slip amplitudes, durations and energy dissipa- 
tion, hinting at the possible presence of scale free distributions. In this paper 
we report on some results from our study of the stick-slip regime of a sheared 
granular medium, when deviations from peridic motion are strong. We describe 
the experimental apparatus in Sec. [21 from which we extract the statistics of 
the reaction force (Sec. [3|), the distribution of which has some notable features 
(Sec.[3|). Section 5 discusses the statistics of velocity fluctuations and compares 
with the stochastic model and Sec. [5] contains a brief summary. 

2 The experiment 

In order to investigate the stick-slip response of a granular medium to an external 
shear rate we devised and realized a circular channel in which a granular medium 
is sheared by an overhead rotating annulus which bears on the medium (Fig.[T|). 
The plate is driven by a motor via a torsion spring, and is free to move vertically. 
As the motor turns the springs loads and the annular plate is initialUy retarded 
by static friction with the granulate. As soon as the force is sufficient to destroy 
the granular structure (i.e. force chains in the medium, see schematic in Fig.[T]), 
the plates starts to slip, unloading the spring, and slows to a stop when the 
friction force overcomes the spring force. Two optical encoders measure the 
angular difference {ujot — 0) between the motor and the plate, and this signal 
is stored in a pc as a function of time. The resolution on the angle is 3 x 10~^ 
rad, and up to 2 x 10~^s in time (though signals are typically smoothed and 
interpolated to 10 ms or 1 mrad, in order to reduce instrumental noise). The bed 
of glass beads is combed before each series of experiments in order to minimize 
ageing effects and also to remove any possible crystal structures present. This 
preparation procedure allows for an acceptable reproducibility of the results. 
More details on the experimental set-up can be found in [2]. 

3 Statistics of the reaction torque 

By measuring the instantaneous angle, several quantities can be derived. For 
instance, a sample of the resulting angular velocity time series of the plate 4>{t) 
and driving motor is given in Fig. [21 The two curves with clearly distinct 
features correspond to the stick slip and to the sliding regimes. 
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Figure 1: The experiment consists of an annular top plate shearing over a 
granular material in a Couette geometry. 
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Figure 2: Shear velocity as a function of time in the stick-shp (sohd) and shding 
(dotted) regimes. The respective driving velocities are shown by the dashed 
lines. 

The reaction torque t of the granular medium can be derived as a functon 
of time: 

T{t)=K{LUDt-(f>it))+I4>{t) (1) 

where k is the spring constant and / the total inertia of the plate and annexes. 
The torque may also be more usefully expressed as a function of the angle: t(^). 
Experiments were performed at various driving velocities and springs constants, 
the variation range of which is illustrated in Tab. [1] The time series exemplified 
in Fig. [2] have been exploited to extract important information on the statistical 
properties of the torque r. In particular we have constructed the Probability 
Density Function (PDF) P(r) and the Auto-Correlation Function (ACF). For 
each of the experiments, we have considered the following values of t: i) at the 
detachment (static), ii) at stationary points of velocity (when = 0), and iii) 
globally. The resulting PDF's are shown in [2] and are virtually identical. Their 
most notable feature, the asymmetry, is clearly evident. Different functional 
forms can be used to reasonably fit the curves (e.g. Gumbel, Lognormal, Gamma 
etc.) but none seem to offer any insight into any possible underlying physical 
process [2]. 

The AGF is calculated on the torque signal as a function of angular displace- 
ment, and in shown in Fig. [3] for a single experiment. The straight line is an 
exponential decay of angular constant (/)c = 1 radian (corresponding to 30,000 
times the system's angular resolution, or roughly 100 particle diameters). This 
observation suggests that a model for the granular friction force should take 
explicitly into account the correlations present in the granular medium. 
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Figure 3: The Autocorrelation Function for the torque exerted by the torsion 
spring. An initial exponential delay with characteristic angle 0c — 1 rad (solid 
line) is followed by a more rapid decay. 

4 Stability and universality 

Non-symmetric distributions appear in many natural phenomena, and are usu- 
ally deemed to be an indication of the presence of correlations. In recent years 
some authors [7] have suggested the idea that in many cases most, if not all, of 
these distributions pertain to a simple universal kind of PDF, which would take 
the place of the Gaussian distribution when the variables are strongly correlated. 

Actually, by virtue of the central limit theorem, the PDF for the sum of a 
large number of independent variable converges to the Gaussian distribution. 
In Ref. [7] and related work, it is proposed that the observed non-symmetric 
distributions may be limit distributions for sums of strongly correlated variables. 

The distribution of the torque shows a very similar behavior, with a pro- 
nounced asymmetry. This is shown in Fig. 4, where the torque distribution is 
the curve denoted by M = 1 (see below). 

In order to test if the torque PDF is compatible, at least in principle, with 
such a limit distribution, we have computed the PDF of the sum variable 

M 

ZM = J2T{cj),)/M 

i=\ 

where (pi = i ■ A(/)o, (i = 1,2,...), and we compare the distribution for several 
values of M. In our data A0o ~ 1 mrad. In figure |4] we report the distributions 
for M = 1, 16, 256, and 512 (i.e. corresponding to 1 ^ 512 mrad). The 
distributions are rescaled in order to have zero mean and unit variance. 
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Rescaled torque 

Figure 4: Comparison of the sum variable Za/'s PDF with rescaled gamma 
distribution with parameter r = 10. The PDF of the torque corresponds to the 
M=l curve. It presents a pronounced asymmetry. 



Recalling the analysis of the torque distribution, if we assume that the vari- 
ables Ti = T{4>i) are Gamma distributed, i.e. according to: 



1 



T[-^exp(-r,/0) 



then for any i the rescaled variable z ~ t-<t> gj^ould be distributed as 



P{z) = —-- cyiY>{-\/rz - r){^/rz + rY 

r(r) 

Fitting our measures for the torque distribution with a Gamma distribution, 
we find that r « 10. Note also that the variables are not independent: their 
correlation coefficient decays exponentially as: p'*^-'!. A comparison with the 
measured autocorrelation function gives 

p = exp(-A0o/0c) = exp(-l/1000) « 0.999 

Accordingly with previous mathematical studies [51 [S] , the distribution of 
Zm can itself be approximated with a gamma distribution with parameter: 



r(M) 



M 



2p 

1-p 



M 



1-p 



(2) 



(the scale parameter is unnecessary if we compare the rescaled variables). Con- 
sidering our value of p, the variation of tm from A/ = 1 to 512 is between 
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Figure 5: The event duration distribution (inset) shows a power-law superim- 
posed with a peak, here fitted by a WeibuU. The mean and mode of the WeibuU 
peak increase as the square-root of the moments of inertia, while the power-law 
exponent remains constant. 



ri = 10 and r5i2 « 11. In Fig. U] we compare the measured rescaled ditribution 
of Zm with the rescaled gamma distribution (r = 10). 

Figure |4] shows how the rescaled distributions do not change, and this phe- 
nomenon can be interpreted as a form of "weak stability" of the torque distri- 
bution [lOj . A similar mechanism could also play a role in other situations and 
may explain the universal behaviour discussed in Ref. [3 [TT] . 



5 The equation of motion 

The work by Dalton and Corcoran [31 [S] , performed at constant volume, revealed 
the presence of self similar PDFs for the duration and extension of the slip 
events, and for the velocity. Our experiment confirmed the presence of such 
distributions, although in this case the medium can dilate and a peak at large 
duration appears [12] (inset to Fig. [5]); this peak was ascribed to the inertial 
properties of the plate. Figure O shows the peak position as a function of the 
moments of inertia of the top plate. It is clear that the peak displaces as VT 
which is to be expected considering the top plate as a harmonic oscillator with 
frequency T — In the same article [12], we also presented a stochastic 

model quantitatively describing the dynamics of the stick-slip motion. 

In order to find an effective description of the dynamics, an extensive analysis 
of the reaction torque has been carried out. Figure [6] shows the reaction torque 
T as a function of the velocity for a set of slip events, apparently showing a 
completely irregular behaviour. If an average is performed however, then it 
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Table 1: Range of values over which the spring constant, k, and motor angular 
velocity, uj can be varied; in parenthesis values typically used in experiments. 



K{Nm/rad) 


uj{rad/s) 


0.12 (0.36) 


10-5 (0.005) 


0.99 


60 (1) 



is possible to identify a statistically regular behaviour. Not surprisingly it is 
characterized by the usual dynamic instability which is common to solid-on- 
solid friction and which is at the origin of stick-slip motion. An adequate fit for 
the curve is (Fig. [6]): 

Td((/)) = TO + 7(0 - 2^0 ln(l + <j>/uJo)). (3) 

What is important is that it contains just two parameters, 7 the high-velocity 
viscosity and ujo, the location of minimum friction, both of which show little 
dependence on the experiment parameters to, I and n. 

Once the deterministic part of the shear stress has been taken into account 
by the friction law, depending only on the angular velocity, one remains with 
only the stochastic part of the torque (which generally behaves as shown in 
Fig [7]). Analysis of the correlation and of the spectral properties of suggests 
that, with good approximation, it performs a bounded random walk in (j) [12j : 

^ ^ D7j{<j>) - ar,, (4) 

r] being a unit white noise and D and a measured from the power spectrum 
which should have the form 

With these few assumptions (Eqs. ([3]), (P), Q) we are in a position to 
simulate the motion of the top plate. The results have been outlined in T2] and 
the model shows excellent agreement with experimental data such as the slip 
durations and sizes, and the velocity distributions. 

Notably, it is quite interesting that an almost identical model as that de- 
scribed here has been used to explain the origin of Barkhausen noise in magnetic 
hysteresis cycles jl3j. where magnetic domain walls respond intermittently to 
a varying external field. Furthermore, through the use of this simple model, 
chaotic granular stick-slip systems, heretofore notoriously intractable, have been 
opened to analysis and characterisation. 



6 Summary and conclusions 

The stick slip regime of a granular medium subject to shear can be statisti- 
cally well characterized in terms of the probablity distribution of stress (torque), 
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Figure 6: The torque as a function of the instantaneous velocity for each shp 
event, their average (circles) and the fit to Eq[3] (thick solid line). 




Angular displacement (rad) 

Figure 7: The random component of the torque after the deterministic part, 
dependent only on the instantaneous velocity, has been removed. The upper 
curve shows the velocity as a function of the angle. Events are delineated by 
moments when the velocity reaches zero. 
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velocities, slips durations and extensions, etc. 

We measure the torque distribution, which is a physically relevant observ- 
able, and we observe non Gaussian, asymmetric fluctuations. This distribution 
is similar to those encountered in many different phenomena involving corre- 
lated variables and has been proposed as a sort of 'universal' distribution for 
sums of correlated variables. We show that our results seem to exhibit some 
degree of stability of the sum variable. However, this result can be explained by 
simply hypothesizing Gamma-distributed torque values, and so a true universal 
limit distribution is not implied, particularly given the absence of a favourable 
convincing theoretical argument. 

In our work we show that, despite the complexity of the mesoscopic inter- 
actions between the grains, the statistics of fluctuations is quantitatively repro- 
duced by a minimal model, based essentially on the Brownian nature of the 
frictional granular forces, accounting for the spatial correlation of the friction. 

Both the experiment and the model exhibit broad slip duration distributions 
which feature a characteristic peak. In this work we present an analysis of 
this peak and show, as expected, that its position is specified by the resonant 
frequency of the plate and spring. 

This work was supported by the Italian FIRB project RBAZ01883Z and the 
EU project "TRIGS", c.n. 043386, COFIN 2005. 

References 

[I] R. A. Bagnold, Proc. R. Soc. A 295, (1966) 219-232. 

[2] F. Dalton, F. Farrelly, A. Petri, L. Pietronero, L. PitoUi and G. Pontuale, 
Phys. Rev. Lett. 95, (2005) 138001-1-4. 

[3] F. Dalton and D. Corcoran, Phys Rev. E 63, (2001) 61312-1-4 

[4] S. Nasuno, A. Kudrolh, A. Bak and J. P. Gollub, Phys. Rev. E 58, (1998) 
2161-2171. frictional granular forces, accounting for the spatial correlation of 
the friction. 

[5] P. A. Thompson and G. Grest, Phys. Rev. Lett. 67, (1991) 1751-1754. 

[6] F. Dalton and D. Corcoran, Phys Rev. E 65, (2002) 31310-1-5 

[7] See e.g. Europhys. Lett. 57 (2002) 310-314 and references therein 

[8] S. Kotz, J. Nemnann, J.Gcophys.Rcs. 68, 3635-3641 (1963). 

[9] S. Kotz and J. W. Adams, The Annals of Mathematical Statistics, 35, 277- 
283, 1964. 

[10] A. Renyi, Calcul des probabilites (Dunod, Paris 1966) 327-331. 

[II] I. M. Janosi and J. A. C. Gallas, Physica A 271 (1999) 448-457. We are 
grateful to the authors for sharing with us their data. 



10 



[12] A. Baldassarri, F. Dalton, A. Petri, S. Zapped, G.Pontuale and L. 
Pietronero, Phys. Rev. Lett. 96, 118002 (2006) 118002-5. 

[13] S. Zapperi, P. Cizeau, G. Durin, and H. E. Stanley, Phys. Rev. B 58, 6353 
(1998); B. Alessandro, C. Beatrice, G. Bertotti, and Ate. Montorsi, J. Appl. 
Phys. 68, 2901 (1990). 



11 



